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Following my previous proposal that two-particle correlation functions can be 
used to resolve the minijet contribution to particle production in minimum biased 
events of high energy hadronic interactions, I study the pj and energy dependence of 
the correlation. Using HIJING Monte Carlo model, it is found that the correlation 
c{4>i,4'2) in azimuthal angle (j) between two particles with pjj > p™* resembles much 
like two back-to-back jets as p™* increases at high colliding energies due to minijet 
production. It is shown that c(0, 0) — c(0, vr), which is related to the relative fraction of 
particles from minijets, increases with energy. The background of the correlation for 
fixed p™* also grows with energy due to the increase of multiple minijet production. 
Application of this analysis to the study of jet quenching in ultrarelativistic heavy 
ion collisions is also discussed. 
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I. INTRODUCTION 

Minijets with px ^ few GeV/c are commonly believed to become increasingly 
important in liadronic interactions for energies beyond the CERN Intersecting Storage 
Ring (ISR) energy range. It has been suggested by many authors that they are 
responsible for not only the global properties, like the rapid increase of the total cross 
sections[l]-[7] and the average charged multiplicity[7, 10], but also the local correlation 
and fluctuations[6]-[10] of multiparticle production in high energy and pp collisions. 
They have also been estimated to play an important role in ultrarelativistic heavy 
ion collisions [11]. 

Another particularly interesting feature of minijets is that they seem to give a 
natural explaination[12] to the apparent "flow" effect in high energy pp collisions. 
This "flow" effect is the observation in experiments[13] at Fermilab Tevatron Collider 
energy that the average transverse momentum of charged particles increases with the 
total multiplicity of the events and the increase is stronger for heavy particles than 
the light ones. However, this observation is surprisingly in coincident with the results 
of an equilibrated quark gluon plasma (QGP)[14] in which the common transverse 
flow velocity from the collective expansion gives heavy particles larger transverse 
momentum than the pions. The speculation was made even more plausible by Levai 
and Miiller's finding[15] that there is no time for the baryons to equilibrate with the 
pions during the expansion of a hadronic fireball. Other models such as the string 
fusion model[16] can also explain the observed phenomenon. In order to differentiate 
the minijet picture from the other scenarios, one must find a way to confirm the 
infiuence of minijets on particle production. 

Although minijets are anticipated from the prediction of perturbative QCD 
(PQCD) and many efforts have been made to investigate their phenomenological 
consequences, there exists little direct experimental evidence of their presence in 
hadronic collisions. The only experimental attempt to find minijets is by UAl ex- 
periments [17] at CERN where hadronic clusters with transverse energy Et ^ 5 GeV 
have been identified as minijets and the cross sections are found to be consistent with 
the PQCD prediction. However, this kind of cluster-finding method even for inter- 
mediate Et values is complicated by the background of random fiuctuations[18]. Eor 
smaller Et^ minijet clusters are overwhelmed by the background fiuctuation. There- 
fore, because of their small transverse momenta, minijets with relatively small px ^ 
2 GeV/c can never be resolved as distinct jets from the large soft background in 
minimum biased events. To avoid the experimental difficulties of reconstructing jets 
with the huge background in high energy AA collisions, it has been suggested[19] that 
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single inclusive spectra of produced particles can be used to study the properties of 
jet production. Similarly, I have proposed in a previous paper[20] that two particle 
correlation functions is useful to study the content of minijets in minimum biased 
events of hadronic interactions. In this paper, I report the detailed study of the px 
and energy dependence of the two-particle correlation in the azimuthal angle (f> and 
the influence of finite rapidity cuts. 

Since particles from jet fragmentation are correlated in both directions of the 
minijets, they have been found[2I] to dominate both the short and forward-backward 
two-particle correlations in pseudorapidity. For two-particle correlation in the az- 
imuthal angle ^, contribution from back-to-back minijets should be strongly peaked 
at both forward (A^ = 0) and backward (A^ = tt) directions. If we calculate the same 
correlation, but for some selected particles whose transverse momenta are larger than 
a certain px cut, the two peaks should be more prominent because these particles are 
more likely to come from minijets. The energy and px dependence of such correlations 
of charged particles must depend on the relative fraction of minijet events and mini- 
jet multiplicity, thus shedding some light on minijet content in these events. On the 
contrary, particles from soft production of fused strings or an expanding quark gluon 
plasma are isotropical in the transverse plane and would only have some nominal 
correlation in the backward direction due to momentum conservation. 

The experimental study of the two-particle correlations is important to clarify 
the controversial issue whether the observed fiavor dependence of the correlation 
between the multiplicity and average transverse momentum at Fermilab Tevatron 
energy [13] is due to minijet production[12], or string interaction[16], or the formation 
of a quark gluon plasma[14]. On the other hand, it can also provide constraints on 
theoretical models as how to combine PQCD hard scatterings with nonperturbative 
soft interactions, how to parametrize the soft processes, and what is the appropriate 
fragmentation scheme of multiple minijets. In the light of the study on intermittency 
in particle production[22], it also helps us to have a better understanding of the 
intermittent fiuctuations in hadronic collisions at collider energies as particles from 
jet fragmentation indeed have been shown to have strong intermittent behavior [23]. 

This study is based on the Monte Carlo model, HIJING[24], which combines a 
simple two-string phenomenology for low px processes together with PQCD for high 
px processes. Since HIJING is a model with a very specific scheme of minijet frag- 
mentation, the results we get are only a qualitative estimate of the effects of minijets 
on the two-particle correlation functions. Though the correlations in pseudorapidity 
T] and azimuthal angle (f> are both important, we limit our discussions only to the 
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latter case. The remaining of this paper is organized as follows. In Sec. II we briefly 
review the HIJING model as discussed in detailed in Ref. [24]. In Sec. Ill we study 
the pt dependence of the two-particle correlations in (f> and their resemblance to large 
Pt jet profiles. Special emphasis is given to the effects of minijets on the back-to-back 
differences of the correlation, c(0,0) — c(0,7r). In Sec. IV, the energy dependence of 
the background of the correlation and its implication on multiple minijet production 
are discussed. The effect of finite rapidity cut is also discussed. Finally, Sec. V 
concludes with a summary and remarks. 

II. THE HIJING MODEL 

The HIJING model[24] has been developed mainly for multiple jet and particle 
production in pp^ pA, and AA collisions at energies ^/s ^ 5 AGeV. The formulation 
of HIJING was guided by the successful implementation of PQCD in the PYTHIA 
model[25] and the need to incorporate a consistent model of soft processes. It thus 
provides a link between the dominant nonperturbative fragmentation physics at inter- 
mediate energies and the perturbative physics at higher collider energies. A detailed 
description of the HIJING model can be found in Ref. [24]. It includes multiple 
jet production with initial and final state radiation along the lines of the PYTHIA 
model[25] and soft beam jets are modeled by quark-diquark strings with gluon kinks 
along the lines of the DPM[26] and FRITIOF[27] models. For nucleus induced reac- 
tions, nuclear shadowing of parton structure functions and final state interaction of 
produced jets are also considered[I9]. 

In the eikonal formalism, the inelastic cross section of pp or pp collisions is given 

by 

= j J25[l _ e-{-soft+^,et)TN{h,s)^^ 

where T]\j{h, s) is the partonic overlap function between two nucleons at impact pa- 
rameter 6, crjet(po, s) is the total inclusive jet cross section calculated from PQCD with 
Pt > Po, and (Tsoft{s) is the corresponding phenomenological inclusive cross section 
of soft interactions. The cross sections for no and j > f number of hard or semihard 
parton scatterings are, 

ao = y" J'6[I - e-"-^'^~(^'^)]e-"^^'^~(^'^) , (2) 

(3) 

with their sum giving rise to the total inelastic cross section, (7^. Following the above 
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formulas, the cross section with at least one hard parton scattering, defined as (Thard^ 
is then 

ahard = I d^il - e-^^«'^~(^'^)]e-^-/'^~(^'^), (4) 
and the average number of parton scatterings per inelastic event is 

{rijet) = crjet{po,s)/a,n- (5) 

Shown in Fig. 1 are the energy dependence of {rijet) and the fractional cross 
section of hard processes, (Thardlc^in for po = 2 GeV/c and (Tsoft = 57 mb. We see 
that as energy increases inelastic events contain more and more hard processes and 
finally are dominated by minijet production. The probabilities for multiple minijet 
production also become prominent at energies above 1 TeV as {rijet) ~ 1- 

Once again, we emphasize here that po and (Tsoft in HIJING model are phe- 
nomenological and model dependent parameters separating PQCD at high px from 
the nonperturbative low px regime. Our guideline for choosing the value of po is that 
it should be much larger than Aqqj^ ~ 200 MeV/c in order to apply PQCD but 
low enough to permit the simplest possible model for the nonperturbative dynamics. 
Choosing Po = 2 GeV/c allows us to use a constant (Tsoft = 57 mb to reproduce the 
observed energy dependence of the total, elastic and inelastic cross sections in pp 
and pp collisions[10, 21]. With these parameters, it has been shown[21] that a simple 
two-string model suffices to account for the soft dynamics. A complete test of the 
HIJING model has been reported in Ref. [21]. The model has been successful to 
provide a consistent explanation of not only the energy dependence of pseudorapidity 
distributions, moderate inclusive spectra, and the violation of Koba-Nielsen-Olesen 
(KNO) scaling, but also the fiavor and multiplicity dependence of the average trans- 
verse momentum[12] in pp and pp collisions in a wide energy range -^s = 50-1800 
GeV. In addition, it is also consistent with the available data on pA and AA collisions 
at moderate energies -^s < 20 AGeV[24]. 

III. TWO-PARTICLE CORRELATION 

Experimental study of particles from the fragmentation of high px jets has shown 
that particles with larger average transverse momenta are much concentrated in both 
directions of the back-to-back jets. The widths of these jet profiles are about 1 in both 
pseudorapidity rj and azimuthal angle (f> as has been determined from the calorimetric 
study of high px jets. It is apparent that similar trends should also hold for particles 
from minijets and they must infiuence the two-particle correlation functions. It has 
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already been shown in Ref. [21] that minijets are the dominant mechanism underlying 
both the enhanced short range and forward-backward correlations in rj. In this paper 
I demonstrate that the energy and px dependence of two-particle correlations in the 
azimuthal angle (f> can be used to study minijets in minimum biased events. One 
can expect that with higher p™*, the correlation must look similar to the high px jet 
profiles. Moreover, the energy dependence of the correlation for fixed p™* should be 
related to the increase of the fraction of minijet events and the average number of 
minijets as shown in Fig. 1. 

The normalized two-particle correlation functions are defined as, 



where p[(f>) is the averaged particle density in (f> and p[(f>i^(f>2) is the two-particle 
density which is proportional to the probability of joint particle production at (pi 
and (f>2- Both p[(f>) and p[(f>i^(f>2) are integrated over the whole rapidity range. The 
limitation of a restricted rapidity range will be discussed in the next section. In 
actual calculations, one direction in the transverse plane is randomly selected as 
^1 = for all the events and then both the averaged single inclusive density p{A(f>) 
and the two-particle joint density /j(0, A^) are calculated. Because the collisions 
are symmetrical in the transverse plane perpendicular to the beam direction, p{A(f>) 
should be a constant equal to the total averaged multiplicity of the selected particles 
divided by 27r. Similarly, /j(0, A^) is also independent of the direction which is chosen 
as ^1 = 0. In fact, one can even choose the direction of each particle as = to 
calculate c(0,A^) and then average it over all particles in each event. In this case 
one can get better statistics for a limited number of events. Since I have the luxury 
to produce many events through the Monte Carlo simulation, I instead try to average 
over as many events as we can to achieve the best statistics possible over global 
quantities such as the number of jets per event with relative large px- 

Shown in Figs. 2-5 are calculated results on two-particle correlation functions 
in pp and pp collisions at different energies. To study the px dependence, the cor- 
relations are obtained for selected particles with different transverse momentum cut 
px > p™*=0 (long-dash-short-dashed histograms), 0.5 (dotted histograms), 1.0 (solid 
histograms), 1.5 (short-dashed histograms), and 2.0 (dot-dashed histograms). At low 
energy -^s = 50 GeV (Fig. 2) where minijet production is not very important as seen 
from Fig. 1, there is very little correlation between two particles with A^ < 27r/3 ex- 
cept some small increase at A^ ~ for the case of px > 1.0 GeV/c. The correlation 
at A^ ~ TT is totally due to the momentum conservation which increases with larger 




(6) 
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Pt cuts. As was pointed out in Ref. [21], minijet contribution to particle production 
at the highest CERN ISR energies -^s 50 GeV is just beginning to emerge. This is 
again illustrated here by the slight increase of the correlation at ~ for particles 
with Pt > 1.0 GeV/c. 

As we can see from Fig. 1, minijet production increases very fast with energy 
and becomes a significant part of the hadronic interactions above -^s ^ 100 GeV. 
Therefore, the two-particle correlation functions in Fig. 3 for -^s = 200 GeV are very 
different from those at -^s = 50 GeV in Fig. 2. We observe that there is apparently 
strong correlation between two particles with px > 0.5 GeV/c at both ~ and tt 
forming a valley at A^ ~ 7r/3. This feature is very similar to the high px jet profiles as 
functions of (f> relative to the triggered jet axis[28]. However, due to the relative large 
fraction of particles from soft interactions where momentum conservation dictates 
the two-particle correlation at A^ ~ tt, the resultant c(0,A^) at A^ ~ tt is still 
larger than at A^ ~ 0. Whereas in high px jet profiles, the peak of particle fiow 
toward the triggered jet (A^ = 0) is always narrower but the height is about the 
same as in the backward direction(A^ = tt). The dominance of particles from soft 
interactions at this energy is clearly demonstrated in the case of p™* = where 
c(0, A^) monotonically increases from 0.2 at A^ = to 0.3 at A^ = tt. 

As we increase the energy to 1.8 and 6 TeV (Figs. 4, 5), the correlation functions 
with no Pt cut become fiatter over the whole angular range due to the increasing 
fraction of particles from minijets relative to the soft processes. The correlation 
functions with higher px cuts look more and more similar to the high px jet profiles 
with c(0, 0) ~ c(0, tt) and narrower peaks at A^ ~ 0. We know from Fig. 2 that 
particles from soft production have some strong correlation only at A^ ~ tt due 
to momentum conservation while particles from minijets, on the other hand, have 
enhanced correlation both at A^ ~ and tt. The actual shape of the correlation 
function then depends on the relative fractions of particles from minijets and soft 
processes. Looking through Figs. 2-5, we can clearly see that the height of the 
correlation functions at A^ ~ tt gradually becomes equal to that at A^ ~ as 
shown in Fig. 6 by the increase of the difference c(0, 0) — c(0, tt) with energy. It 
approaches zero for any fixed p™* at higher energies. This energy dependence of 
c(0, 0) — c(0, tt) is a direct consequence of the increase of minijet production with 
energy in pp and pp collisions. Thus, it has been demonstrated that the study of the 
energy and px dependence of two-particle correlation function c(0, A^) can provide 
us with information about the minijet content in minimum biased events. 
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IV. EFFECTS OF MULTIPLE 

MINUETS AND FINITE RAPIDITY CUT 

When the average number of produced minijets is larger than one (Fig. 1) at high 
coUiding energies, multiple minijet production becomes important [10, 24]. Since these 
minijets are independently produced, they should increase the background to the two- 
particle correlation. The very extreme cases are central events of high energy heavy 
ion collisions where hundreds of minijets are produced. The two-particle correlation 
will be completely flat when the valleys we see in Figs. 2-5 for pp and pp collisions 
are filled up by the increased background. In hadronic collisions this rarely happens 
and the increase of the background can be considered as an implication of multiple 
minijet production. 

Plotted in Fig. 7 is the difference c(0,0) — c(0,7r/3) between the correlation at 

= and the background in the valley = 7r/3 for each fixed p™* as a function 
of y^. It grows first with energy due to the increase of minijet events and the relative 
small number of minijets per event. However, as the probabilities of multiple minijet 
production become increasingly important, the difference then decreases with energy. 
The larger the p™*, the sooner the decrease begins. 

Unlike high px back-to-back jets which are both kinematically bounded to the 
central rapidity region, a pair of minijets can be easily produced with large rapidity 
gap between them. When we trigger one minijet in a limited rapidity window, the 
other one which is produced in the same parton scattering often falls outside the 
rapidity window. Therefore, if we calculate two-particle correlation for particles in 
a limited rapidity range, minijet contribution to the backward correlation (A^ = tt) 
will mostly drop out while the contribution to forward correlation (A^ = 0) still 
remains. Indeed, as shown in Fig. 8, the forward correlation at A^ = for particles 
in \r]\ < 1 is very strong, but the backward correlation at A^ = tt is drastically 
reduced as compared to the correlation pattern in the full rapidity range in Fig. 2-5. 
Furthermore, due to strong short range two-particle correlation in rapidity [21], the 
forward correlation at A^ = is enhanced by restricting particles to < 1. At 

=1.8 TeV, I find that the enhancement of backward correlation at A^ = tt due 
to minijets becomes important only when the rapidity window is \r]\ ^ 2. 

V. SUMMARY AND REMARKS 

It is demonstrated in detail in this paper that two-particle correlation functions 
in the azimuthal angle (f> are useful for the study of minijets and their contribution to 
particle production in the minimum biased events of hadronic interactions. Because 
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particles from jet fragmentation are more concentrated along the jet axes, they result 
in enhanced two-particle correlation in both the triggered (A^ = 0) and the backward 
(A^ = tt) direction. On the other hand, particles from isotropic soft production only 
have some nominal correlation in the backward direction due to momentum conser- 
vation. Therefore, the energy dependence of the relative height c(0, 0) — c(0, tt) of the 
correlation in the two opposite directions is directly related to the minijet content in 
the particle production and increases with the colliding energy. Furthermore, mul- 
tiple minijet production becomes important at high energies and they increase the 
background relative to the enhanced correlation at A^ = and tt. 

It is important to keep in mind that the results we obtained via HIJING are 
only intended to show qualitatively the effects of multiple minijets on the energy 
and pt dependence of two-particle correlation in (p. The quantitative magnitudes of 
the correlations must depend on the specific scheme of minijet fragmentation and 
the color fiow of multiple minijets[24] in HIJING. It is important to study in the 
future the sensativity of two-particle correlation functions to the schemes of minijet 
fragmentation and the color fiow structure of multiple minijets. The final schemes 
which survive the tests on single inclusive spectra have also to be constrained by 
the experimental measurements of the energy and px dependence of two-particle 
correlation functions. 

I have emphasized the implications of the energy and px dependence of the two- 
particle correlation functions on minijet production. For experiments at a fixed en- 
ergy, one can also investigate the contribution of minijets to the two-particle correla- 
tion function by comparing the direct measurements to the calculation from the same 
data sample but with randomized azimuthal angle (f> for each particle. Due to this 
randomization of ^, the correlation calculated from the manipulated data should be 
much fiatter than the real one for any given p™*. It is also important and interesting 
to study the px dependence of two-particle correlation functions in pseudorapidity rj 
and the correlation relative to a triggered particle which has the highest px in each 
event [20]. 

As discussed in Sec. IV, the two-particle correlation in (f> in ultrarelativistic heavy 
ion collisions should be much fiatter than in pp and pp collisions for any fixed p™* 
due to the huge number of minijet production. However, if we increase p™* to a much 
larger value (e.g. 4-6 GeV/c), the background to the correlation becomes smaller 
due to small cross sections of large px jets. Since the single inclusive spectrum at 
large and moderate px has been shown to be sensative to jet quenching[I9], two- 
particle correlation is apparently more relevant to the study of jet quenching and the 
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implication on energy loss mechanism it conveys. 
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FIGURES 



FIG. 1. Fractional cross section cfhardl<^in for at least one hard scattering and average 
number of minijet production {rijet) per event as functions of ^ys. 

FIG. 2. The correlation functions c(0, A(^) between two particles with fixed px cutoff 
vs Acj) in pp collisions at ^/s = 50 GeV. 

FIG. 3. The same as Fig. 2, except for pp collisions at ^/s = 200 GeV. 

FIG. 4. The same as Fig. 2, except for pp collisions at ^ys = 1.8 TeV. 

FIG. 5. The same as Fig. 2, except for ^ys = 6 TeV. 

FIG. 6. The energy dependence of the difference c(0, 0) — c(0, vr) between the correla- 
tion at Acj) = and vr for each fixed px cutoff. 

FIG. 7. The energy dependence of the difference c(0,0) — c(0,7r/3) between the cor- 
relation at Acj) = and the background at vr/S for each fixed px cutoff. 

FIG. 8. The same as Fig. 4, except for charged particles in a limited rapidity range of 
kl < 1- 
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